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The QCD free energy can be studied by dimensional reduction to a three-dimensional (3d) ef- 
fective theory, whereby non-perturbative lattice simulations become less demanding. To connect 
to the original QCD a perturbative matching computation is required, which is conventionally 
carried out in dimensional regularization. Therefore the 3d lattice results need to be converted to 
this regularization scheme as well. The conversion must be carried up to 4-loop order, where the 
free energy displays an infrared (IR) singularity. We therefore need a regulator which can be im- 
plemented both on the lattice and in the continuum computation. We introduce a mass regulator 
to perform Numerical Stochastic Perturbation Theory computations. Covariant gauge is fixed in 
the Faddeev-Popov scheme without introducing any ghost fields. 
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1. Introduction 

As is well known, finite-temperature QCD seems to show two different phases: it is confining 
at low temperatures (the realm of mesons and baryons) while asymptotic freedom and a quark- 
gluon plasma are expected to appear in the high-temperature regime. A good observable to witness 
the change is the QCD free energy density, given essentially by the familiar Stefan-Boltzmann law 
of blackbody radiation, multiplied by the number of light effective degrees of freedom. 

To study the free energy density requires different methods in different regimes. At low tem- 
peratures the problem has to be treated with numerical lattice simulations, while at high tempera- 
tures perturbation theory should allow at least for some progress, given that the coupling constant 
g is small. Nevertheless, even for small g, certain coefficients in the weak-coupling expansion do 
remain non-perturbative [j|], and can only be determined with numerical techniques. 

In the high-temperature regime, the theory contains three different momentum scales [fj], 
namely T (hard modes), gT (soft modes) and g 2 T (ultrasoft modes). The contribution of each 
of these modes is best isolated in an effective theory setup. This is accomplished via dimensional 
reduction [Q, by integrating out the hard and soft modes to obtain a 3d pure Yang-Mills SU(3) 
theory ("MQCD"). MQCD can then be analysed on the lattice and the results can be added to the 
various perturbative contributions to obtain the complete answer. 

To add the MQCD lattice results to the perturbative ones, we need to change regularization 
scheme from lattice to dimensional regularization. To this aim, a matching between lattice and 
continuum computations is needed and this is achieved by means of Lattice Perturbation Theory 
applied to MQCD. The strategy we adopt for this purpose here is the one of Numerical Stochastic 
Perturbation Theory (NSPT) developed in recent years by the Parma group. 



2. The NSPT method 

NSPT relies on Stochastic Quantization [j|] which is characterized by the introduction of an 
extra coordinate, a stochastic time t, together with an evolution equation called the Langevin equa- 
tion, 

where tj (x, t) is a Gaussian noise which effectively generates the quantum fluctuations of the theory. 

The average over this noise is such that, together with the appropriate limit in t, the desired 
Feynman-Gibbs functional integration is reproduced: 

(0[^(x,t)]) n ^ y [D<l>]0[<i>(x)}e- Sm] ■ (2-2) 
For SU(3) Yang-Mills theory, the Langevin equation becomes 



d t U 1 , = -i(VS[Ur,] + ri)Ur 1 , (2.3) 



guaranteeing the proper evolution of variables within the group. 

In this framework, perturbation theory comes into play by means of the expansion 

k 



tf,,(x,0— ►E*o^* ) (*>0, ( 2 - 4 ) 
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where go is the bare gauge coupling. This results in a system of coupled equations that can be 
numerically solved via a discretization of the stochastic time t = nX, where z is a time step. In 
practice, we let the system evolve according to the Langevin equation for different values of T, 
average over each thermalized signal (this is the meaning of the above-mentioned limit ?—►«>), and 
then extrapolate in order to get the % = value of the desired observable. This procedure is then 
repeated for different values of the various parameters appearing in the action. 



3. Mass as an IR regulator 

As stated above, the quantity we are interested in is the contribution to the QCD free energy 
density / coming from the 3d pure SU(3) theory. On the lattice, this observable is related to the 
trace of the plaquette (1 — Hp), where Hp = N~ l ReTrP and P is the elementary plaquette, via 



2a d d If 
d(d-l)dpo\T 



[\-U p ) = — — — — | - ) , (3.1) 



with bare lattice coupling j8q = 2N C / (a go)- The outcome can be expanded in powers of /3o as 



(1 - n " > ^ + l + l + l +0(fir5) - 



The determinations of the first three coefficients in the present setting have been discussed in 
Ref. |0| . The non-perturbative value of the whole quantity has been determined with lattice si- 
mulations in Ref. Terms of 0(j8 o ~ 5 ) disappear in the continuum limit, thanks to the super- 
renormalizability of the theory. Thus only the fourth order coefficient is missing at the moment. 

As shown parametrically in Ref. [jl]] and explicitly in Refs. [Qj, the coefficient q is actually IR 
divergent, and consequently an appropriate regulator must be introduced for its determination. In 
a non-perturbative setting this is provided by confinement, while in fixed-order computations one 
could employ a finite volume (as in Ref. [Q]) or a mass. Since the use of a mass is more convenient 
in continuum computations involving dimensional regularization, we need to implement it in lattice 
perturbation theory as well. 

Apart from introducing a mass, we also fix the gauge in order to match the setting of the 
continuum computations. Consequently, the functional integral is given by 

Z = |[D0]Det(-£^D M [0]+m 2 )exp(-5 vv -5 Gf ) = J[D(j)] exp ( - S w - S GP - S PP ) , (3.3) 
where we assume the use of lattice units (i.e. a = 1), and 

S w = /3 £(l-n P ) + §^ £ , (3.4) 

P ^ Vc x,)l,A 



/J 



Tr 



ln(-£^[0]+m : 
n 



(3.5) 
(3.6) 
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where we have followed the conventions of Ref. [|1C|], writing in particular = exp(/0^), 0^ = 
<pf^T A , with the normalization Tr[T A T B ] = 8 AB /2. Moreover m is the common gluon and ghost 
mass, a is the gauge parameter, and is the discrete Faddeev-Popov operator, given by jlO| ] 



DM 



1 + _d>„ <E> 2 d> 4 <T> b + 0(<£> ) 

2 12 " 720 " 30240 ^ + { ^ 



d* + i*r, (3.7) 



with <J>^ = (j)^F A , where [F A ]gc = —if ABC are the generators of the adjoint representation. 

To treat the Faddeev-Popov determinant as a part of the action means that, because of the 
Langevin equation, one has to face the quantity VS FP = — VTr[lnfi] = — Tr[VBB~ l ] with B[(j>] = 
— L/i^/j^i[0] +m 2 . We perform the inversion as in Ref. [pi]], while the trace is computed by 
means of sources in the usual way. 

The global strategy is then to perform simulations with lattices of different sizes at fixed mass 
in order to extrapolate to infinite volume and, afterwards, to repeat this procedure for other values 
of the mass. At this point, after subtracting the expected logarithmic divergence, one extrapolates 
to zero mass, obtaining the needed fourth order coefficient. It is crucial to take the infinite-volume 
limit before the zero-mass one because, by performing the limits in the opposite order, the final IR 
regulator would be the volume and not the mass as we want. 



4. First (benchmark) results 

So far, the statistics we took are not sufficient to carry out the infinite- volume and zero-mass 
limits for the fourth order coefficient £4, but it is already possible to crosscheck the reliability of 
the general method. As a first test, we compare the 1-loop numerical results for the trace of the 
plaquette for the various masses with the known analytic values. As shown in Fig. 1 for a lattice 
extent L = 5, the agreement between the numerical values and the analytic curve is very good. 

A second check could consist of extrapolating at fixed lattice extent to zero mass, to see if one 
recovers the already known coefficients [^]. Figs. 2-5 show these extrapolations for a lattice extent 
L = 7: the fitting curve is a polynomial in m 2 (the most naive choice) and it seems to approach the 
expected result (the point at m = 0) very well for all the loop orders. The numerical values are 
given in Table 1 . Both of the mentioned checks are well satisfied also for the other lattice extents 
that we have employed so far. 



Loop 


Result from a fit to m = 


Direct measurement at m = 


1 


-2.6594(17) 


-2.6580(8) 


2 


-1.9166(63) 


-1.9095(30) 


3 


-6.304(37) 


-6.307(21) 


4 


-28.43(27) 


-28.68(15) 



Table 1: Comparison of the zero-mass extrapolations with the known results ||7|] (lattice extent = 7). 



As for the 4-loop order, Fig. 6 shows the behavior with respect to the lattice size at fixed mass: 
the result seems to stabilise towards the infinite-volume value in the way one would expect. Once 
a few more lattice sizes are available and similar extrapolations can be carried out for all masses, 
we will finally be in a position to carry out the mass extrapolation that is our ultimate goal. 



189/4 



Four-loop plaquette in 3d with a mass regulator 



Christian Torrero 



-1.6 




0.5 1 1.5 

Mass 



Figure 1: The 1-loop trace of the plaquette vs. 
mass (for L = 5): the numerical results (blue 
dots) agree with the analytical red curve. 
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Figure 3: The 2-loop trace of the plaquette (for 
L = 7), together with a polynomial fit in m 2 . 
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Figure 2: The 1-loop trace of the plaquette vs. 
mass (for L = 7): the fitted curve in red ap- 
proaches the expected value at m = 0. 
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Figure 4: The 3-loop trace of the plaquette (for 
L = 7), together with a polynomial fit in m 2 . 



5. Conclusions and prospects 

It is worth stressing once again that our approach has successfully passed the reliability checks 
we adopted: known zero-mass limits are reproduced through an extrapolation, and the volume 
dependence at a fixed mass appears to disappear once the dimensionless combination mL, where m 
is the mass and L the lattice extent, is large enough. 

In order to obtain the asymptotic large-volume value at a fixed mass, it is still necessary to 
collect more statistics on bigger lattices (for example, L = 12 and 14) at least for the two or three 
smallest masses. Then, the fitting function should be a combination of a negative exponential and 



polynomials in mL, as explained for instance in Ref. ni2Q. 

After subtracting the logarithmic divergence from the fitted infinite-volume values, the sub- 
sequent extrapolation to zero mass does not appear to be troublesome, given that tests with lower 
loop orders have produced good results so far. 
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